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DYNAMICS OF SPACECRAFT CONTROL LABORATORY 
EXPERIMENT (SCOLE) SLEW MANEUVERS 

Y. P. Kakad 

Dept, of Electrical Engineering 
University of North Carolina at Charlotte 
Charlotte, NC 28223 


SUMMARY 

This is the first report of a set of two reports on the dynamics and control of 
slewing maneuvers of NASA Spacecraft Control Laboratory Experiment (SCOLE) 
article. In this report, the dynamics of slewing maneuvers of SCOLE are developed 
in terms of an arbitrary maneuver about any given axis. The set of dynamical 
equations incorporate rigid-body slew maneuver and three-dimensional vibrations 
of the complete assembly comprising the rigid shuttle, the flexible beam, and the 
reflector with an offset mass. The analysis also includes kinematic nonlinearities of 
the entire assembly during the maneuver and the dynamics of the interaction 
between the rigid shuttle and the flexible appendage. The final set of dynamical 
equations obtained for slewing maneuvers are highly nonlinear and coupled in 
terms of the flexible modes and the rigid-body modes. 

The equations are further simplified and evaluated numerically to include the 
first ten flexible modes and the SCOLE data to yield a model for designing control 
systems to perform slew maneuvers. 
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1. INTRODUCTION 

The primary control objective of the Spacecraft Control Laboratory Experi- 
ment (SCOLE) is to direct the RF Line-Of-Sight (LOS) of the antenna-like 
configuration towards a fixed target under the conditions of minimum time and 
limited control authority [1]. This problem of directing the LOS of antenna- like 
configuration involves both the slewing maneuver of the entire assembly and the 
vibration suppression of the flexible antenna-like beam. The study of ordinary 
rigid-body slew maneuvers has received considerable attention in the literature 
[2,3] due to the fact that any arbitrary large-angle slew maneuver involves 
kinematic nonlinearities. This is further complicated in the case of SCOLE by vir- 
tue of a flexible appendage deployed from the rigid space shuttle. The dynamics of 
arbitrary large-angle slew maneuvers of SCOLE model are derived in this report as 
a set of coupled equations with the rigid-body motions including the nonlinear 
kinematics and the vibratory equations of the flexible appendage. 

The dynamical equations of slewing maneuvers of this large flexible spacecraft 
are developed by writing the total kinetic and potential energy expressions for the 
entire system. The energy expressions are further utilized in formulating 
Lagrange’s equations which are expressed in terms of non-generalized co-ordinates 
using an inertial co-ordinate system and a body-fixed co-ordinate system at the 
point of attachment of the flexible beam to the shuttle. The generic model used for 
this analysis consists of a distributed parameter beam with two end masses. The 
three dimensional linear vibration analysis of this free-free beam model with end 
masses [4] is incorporated together with rigid-slewing maneuver dynamics which 
are written in terms of four Euler parameters [5] and angular rotation about an 
arbitrary axis of rotation to yield the final set of highly nonlinear and coupled 
equations. In the derivation of the equations, it is assumed that the vibratory 
analysis is for small motions. 
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2. LIST OF SYMBOLS 

qjz ) Position vector of mass element on the beam from the point 

of attachment 
B Damping matrix 

C Inertial frame to body-fixed frame transformation 

c_ Position vector from the point of attachment to the mass center 

of the beam 

D Mass density of the beam 

d_{z ,t ) Displacement vector of mass element in the body-fixed frame 
E Modulus of Elasticity 

£, (f ) Force applied at the orbiter mass center 
£(f ) Force applied at the reflector mass center 
Ga (f ) Moment applied about the orbiter mass center 
Modulus of rigidity for the beam 

I Beam cross section moment of inertia 

I x Beam cross section moment of inertia, roll bending 

I y Beam cross section moment of inertia, pitch bending 

I I Mass moment of inertia matrix of the shuttle 

1 2 Mass moment of inertia matrix of the reflector 

J Mass moment of inertia matrix of the beam 

L The Length of the beam 

M Angular velocity vector transformation 

to Total mass of the flexible beam 

to i Mass of the orbiter 

to 2 Mass of the reflector 

n The maximum number of modes considered 

q t Generalized coordinates 
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R 

Position vector of the mass center of the orbiter in the inertial 
frame 

r_ 

Position vector f rom the orbiter mass center to the point of 
attachment 

r x 

x co-ordinate of the reflector mass center in the body-fixed 
frame 

r y 

y co-ordinate of the reflector mass center in the body-fixed 
frame 

T 

Total Kinetic Energy 

U 

Total Potential Energy 

u x (z ,t ) 

The beam deflection in x direction referred to the body- 
fixed frame 

U y C Z ,t ) 

The beam deflection in y direction referred to the body- 
fixed frame 

u^(z ,t ) 

The torsional deflection about z axis in the body- 
fixed frame 

V 

Velocity vector of the mass center of the orbiter in the body- 
fixed frame 

| & 

Velocity vector of the point of attachment in the body- 
fixed frame 

P 

Mass per unit length of the flexible beam 

X 

Vector representing the axis rotation during the slew 
maneuver 

<t>xl 

i th Eigenfunction corresponding to u x 

<t>yi 

i th Eigenfunction corresponding to Uy 

4>1>i 

i th Eigenfunction corrsponding to 

0 

The attitude of the orbiter in the inertial frame 


Slew Angle 

63_ 

The angular velocity of the orbiter in the inertial 
frame 

a 

The angular velocity of the reflector in the inertial 
frame 

C 

Damping ratio 


I 

I 
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3. ANALYTICS 

Co-ordinate Systems 

The motion of SCOLE assembly when considered as a rigid body in space has 
six dynamic degrees of freedom: three of these define the location of the mass 
center, and three define the orientation (attitude) of the body. The motion of this 
rigid body is goverened by newtonion laws of motion expressed in terms of 
changes in linear momentum and angular momentum. These relationships are 
valid only when the axes along which the motion is resolved are an inertial frame 
of reference [9,10]. To define the orientation of the orbiter in space, a set of orthog- 
onal axes fixed in the body is utilized. Then the attitude of the orbiter is defined in 
terms of the angles ( 0 1 , 0 2 »^ 3 ) between the body- fixed axes and the inertial co- 
ordinate axes. The body-fixed frame origin is located at the point of attachment of 
the flexible appendage with the rigid shuttle for this analysis (Fig. 1). 

The transformation from the inertial frame to the body-fixed frame is given 
by the matrix, C as developed in figure 2 where if 7 , /, k represent the dexteral 
set of orthogonal unit vectors fixed in the body- fixed frame and 0 t is the rotation 
about 7 , 0 2 is the rotation about ~j and 0 3 is the rotation about k. These rotations 
are carried out successively as shown in figure 1 and the matrix C is given as 


C = 


cos0 3 sin0 3 0 
— sin0 3 cos0 3 0 
0 0 1 


cos0 2 0 — sin0 2 
0 1 0 
sin0 2 0 cos0 2 


1 0 0 
0 cos0j sin0j 
0 — sin0j cos0 ! 


Thus C T is obtained as 


( 1 ) 


C r = 


COS0 2 COS0 3 

sin0 1 sin0 2 cos0 3 +sin0 3 cos0 1 
— cos0 ]Sin0 2 cos0 3 +sin0 3 sin0 1 


— cos0 2 sin0 3 sin0 2 

— sin0 x sin0 2 sin0 3 +cos0 3 cos0 1 — sin0 jcos© 2 
cos0 1 sin0 2 sin0 3 +cos0 3 sin0 1 cos0 1 cos0 2 


12 ) 
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In order to completely define the attitude (orientation), it is needed to relate 
the rotation angles Oj, 0 2 , and 0 3 to the angular velocity components (a> 1 , co 2 , o> 3 ) 
of the orbiter. One way of obtaining the required relations is via body-three angles 
method [5] which was utilized in developing C matrix in equation (1) and these 
relations are 


0i = (<o 1 cos0 3 — o> 2 sin0 3 )/cos0 2 

0 2 = (&> 1 sin0 3 +o> 2 cos0 3 ) (3) 

0 3 = (— O) 1 cos0 3 +c«) 2 sin0 3 )tan0 2 +<i)3 

Thus, the angular velocity of the orbiter can be obtained in the inertial frame 
by means of the following transformation 


co = M t 9 

where the transformation M T is given as 


(4) 


M t = 


cos0 2 cos0 3 sin0 3 0 
— cos0 2 sin0 3 cos0 3 0 
sin0 2 0 1 


(5) 


Although the body-three angles method is used here for obtaining the 


transformations C and M , there are three other methods which can be used to 
obtain the same transformations. A detailed discussion of all the methods is given 
in reference [5] and a summary of the transformations using the remaining three 


methods is given in the Appendix. 


Kinetic Energy 

If the position vector of the mass center of the orbiter in the inertial frame 
(Fig. 3), R, is given as 
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R = 


Rx 

Ry 

Rz 


then the velocity of the mass center in the inertial frame is 


Z 7 0) = 


Rx 

Ry 

Rz 


This velocity can be transformed in the body-fixed frame as 


( 6 ) 


(7) 


zco= C 


Rx 

Ry 

Rz 


The velocity of the point of attachment in the body-fixed frame is 


( 8 ) 


& = V+wxr (9) 

where r. is the vector from orbiter mass center to the point of attachment. 

Defining the position vector (Fig. 4), a_, of a mass element on the beam from 
the point of attachment (origin of the body-fixed frame) before deformation as 


a = 


0 

0 


z 


and the displacement vector of this mass element as 


u x (z ,t ) 

Uy (z ,t ) 

0 


( 10 ) 


d (z ,t ) = 


( 11 ) 
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the position vector after deflection is given as a+d_. The kinetic energy in the 
beam [6] is 


T 1 = (1/2) mVjV^ + (1/2) gS [J ]<o - mVo T [c ]ca+ (1/2) f dj ddm 


+ 


Vj fddm+gi J J addrn +(1/2) 


dl 


“x 

ily‘ 

U yf, 


( 12 ) 


where the vector c_ is from the point of attachment to the mass center of the beam 
and if it is assumed that the beam is a thin rod, then it is given as 


c = 


= (1/m ) f adm = 


0 

0 

-LI 2 


(13) 


and using the skew symmetric form for the vector cross product for any two vec- 
tors c. and w (in the same reference frame) as 

C_X = [C]<0 


0 -c, 


c = 


C v 


0 -c x 
c, 0 


also, the moment of inertia matrix is given as 


(14) 


J = (1/3) pL 3 


1 0 0 
0 1 0 
0 0 0 


(15) 


where p is the mass per unit length of the beam. The last term in the equation 
(12) corresponding to torsional motion is given as 
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= (1/2) f u x ' Uy' 


(1/2) f [ u x ' Uy' U f 
H2(pds)s 2 0 


dl 


u T 


0 

0 


Uy 

U + 

0 

H2(pds)s 2 0 
0 0 


u. 


Uy' 

Urf, 


(16) 


The kinetic energy equation (12) can be simplified as 


Ti = (1/2) pLVjVg +(1/6) pL 3 as?+a>i -pLVjc.to+pL £ q t : 

' J i = i 


where 


+Vjgt+(£>? /3+(l/4) p 


LP5i4i 2+ LP6i4i 7 

i = 1 i = 1 


(17) 


u x = £ <2> x < Cs )& (? ) 

i = l 
n 

Uy = £ 0 yi (S )?j ( t ) 

( = 1 

u x ' - £ 0 xi 'Cs )q t (r ) 

i = l 
n 


Uy' = £ </>y< '($ )?i (f ) 

i = 1 
n 

U^= £ <t>fi (s )?£ (t ) 


i = l 

L 

= f<t>xi(s)ds 

L 


Pli 

Pli- |</>y iCs)<2S 

L 

P3i = fstixi^s^ds 
0 
L 

p4i = ^s4> yi is')ds 

L 

P5i = S^<t> xi ‘^ds 
P6i = J(s<v) 2 ^ 


( 18 ) 
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and 


4(0 = p 


Epi<4 

i=l 

n 

Lp2i?i 

i = i 

o 


4(0 = p 


n 

LPAiti 

i = 1 
n 

]£p3i4 
1 = 1 

0 


(19) 


( 20 ) 


The expressions for p li( p 2 * , p 3i , p 4* , p 5* , and p 6 * are developed as follows. Note 
that 


(t> xi (s ) = A x * sin/3 * s +Z? X * cos/3 * s +C X * sinh/3 * s +D xi cosh/3 * s 


where /3* = 


t* 2 P 

£/ 


1/4 


Since for SCOLE configuration EI X = EI y and = /3* y , £7 and /3* are used 
for both 0 X * (s ) and </> y * (s ). However, this may not be true for other 
configurations. 


Pu = 


ViL 1 


Pli = ftAs^ds 

L 

Pli = f<t>yi(s)ds 

■A xi cos/3* L +B xi sina* +C X * cosh/3* L 
+D xi sinh/3* Z. +A xi —C xi 


(21 A) 


Defining a* = /3*L 
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pu = — I — AriCOSOii+BxiSincXi+CxiCOshai+DxiSinhai+Axi— C xi 

oiiL l 
similarly, 


Pu = 


P 2i = 


— 1 -T \—A yi cosfiiL+Byi sin/3 * A 
/3 iL L l 

+C yi coshj3j L +D yi sinh/Bj L +A yi —C yi 
| — A yi cosoii +B yi sina i +C yi coshar* +D yi sinha* +A yi — C yi 


oil L 


Pa 


PM 


L. 

- fs (}) xi ( s Ids 

L 

= £s(t> yi (s')ds 


and these can be given as 


P 3i = 


A 

sin/3 1 L 

Lcos&iL 

4- R 

cosfiiL ^ L sin/3 t L i 

A xi 

fil 

fit 


(3, 2 + (3, (3? 


K. 

L cosh/3 i L 

sinh/3* A 

4- D 

Asinh/3* A 

cosh/3 t L t i 

u xi 

fii 

&l 

xi 

fit 

S, 2 + e. 2 


Pii 


= A 


xi 


'X i 


A 2 sina f 

L 2 cosaj 

+ B xi 

L 2 cosaj A 2 sina* 

A 2 

a 2 

<*i 

a 2 + 



.L 2 cosha!i 

A^inha* 

4 - T) 

A 2 sinha* A^osha* + L 2 


od \ 

^ U xi 

0<i Oil ^l 


Similarly, 


Pm = 


Ayi 


sin/3 t L 

L cos/3 jL 

4-7? 

cos/3 j L 

Lsin/3jL i 

$1 

0i 


fil 

Ri R? 


'yi 


L cosh/3 jL 

sin/3 jZ, 

+£ ) yi 

Asinh/3j A 

cosh/3* A 

h _l_ 

0i 

0« 2 

0* 

0i 2 

07. 


Pm = 


Ayi 


A 2 sina* 

A 2 cosaj 

4- R 

A 2 cosa e A 2 sinai 

A 2 

<* 2 


-r&yi 

Oil “i 

"«, 2 


'yi 


A 2 coshai 

A 2 sinhar* 

+ D yi 

A 2 sinhai 

L 2 coshoii 

.LL 

°ii 

Oil \ 

Oii 

Oil 



(2 IB) 


(22A) 

(22B) 


(23A) 

(23B) 


(24A) 

(24B) 
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L 

Psi = £ |s0 xi '(s) ds 
P6i = f ) Z ds 


and these can be shown to be 


, 3 ?L 3 i l 


— +I y cos2«,i+(- 


i 2 s, 1 


-)sin2/3,L 


-A, t B xt isin2S,i+ -±—L£l C os20,L-- ±- 

' 

+A xi C xi RiL 2 (cos/3 i L sinh/3 1 L)+ (sin/3 1 L coshj3 1 L ) 

—2 L sin/3 jLsinh^L — (cos/3, L sinh/3 (sin/3 , L cosh/3,L ) 

Pi 

+A xi D xi /3 , Z, 2 (cos/3 , L cosh/3 , Z, )+ (sin/3 , L sinh/3 , L ) 

•2L sin/3 ,L cosh/3 ,L — (cos/3 ,L cosh/3 , L )— (sin/3 t L sinh/3 ,Z- ) +-^- 

Pi Pi 


+BJ ^i-iUcos2/8,L+ 


6 2 2 


2 4/3, 


— B xi C xi /3 i L 2 1 (sin/3 , L sinh/3 , Z, )— (cos/3 , L cosh/3 , L ) | 
+2L cos/3 ,L sinh/3 ,L — -j— (cos/3, L cosh/3, L ) 

) ^ l 

+(sin/3,L sinh/3, Z, ) +-^- 

Pi 

—B xi D xi PiL 2 (sin/3 , Z, cosh/3 , Z, )- (cosj3 , L sinh/3 £ L ) 

+2L cos/3, L cosh/3, Z, — (cos/3 , L sinh/3 , L) 

Pi 

r \ i 

— (sin/3 , L cosh/3 , Z, ) 

+C X ? j 3i^cosh23i^+(-^^+-^)sinh 2 3iZ,--^^- 


Bi-L 2 

2 


-f C X j — “ — cos2j3 * L — -^sinh2j3 * L + — — cosh2j8 * L — ■ 
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+ n 4 

~ ^ XI 


V 6i Ayi 


3, 2 Z 3 


-Ay i B yi 


L 

— -CO! 
2 

L . 
—sir 
2 


/3< Z 2 

, 1 ) 

2 

4/3 £ 3 

r ^ 2 0i 

1 

L 2 

40, 

r 1 _ 

Z 2 /3 £ 

4/3 £ 

2 


£ 2 Z 3 

3 


( 25 ) 


)sin2j3 j Z 
)cos2/3 £ Z 


4/3 j 


+ a r . 


/3,z ; 


(cos/3 j Z sinh/3 £ Z )+ (sin/3 £ L cosh/3 £ Z ) 


— 2L 


sin/3 £ Z sinh/3 £ Z 


1 

/3, 


(cos/3 £ Z sinh/3 £ Z )— (sin/3 £ Z cosh/3 £ Z ) 


+^yi D yi 


0,Z ; 


(cos/3 £ L cosh/3 £ Z )+(sin/3 i Z sinh/3 £ Z ) 


-2Z 


3i 


(cos/3 1 L cosh/3 * Z )— (sin/3 £ Z sinh/3 £ L ) 


sin/3 £ Z cosh/3 £ Z 

1 


0i 


+B yi 2 


(S?L 3 i 


y cos2/3 £ Z + ( )sin2/3 £ Z 


'-®yi Cyi 


PiL 2 


6 2 

(sin/3 £ L sinh/3 £ Z )— (cos/3 £ L cosh/3 £ Z ) 


+2Z 


1 

IS i 


(cos/3 £ Zcosh/3 £ Z )+(sin/3 £ Zsinh/3 £ Z ) 


cos/3 £ Z sinh/3 £ Z 

1 


-B yi D yi 


/3;Z 2 


(sin/3 £ Z cosh/3 £ Z )— (cos/3 £ Z sinh/3 £ Z ) 


+2Z 


J3| 

cos/3 £ Zcosh/3 £ Z 


/Si 


c 2 


(cos/3 £ Z sinh/3 £ Z )— (sin/3 £ Z cosh/3 £ Z ) 
2 


2r 3 


Cyj D yi 


/3 £ Z cosh2/3 £ Z +( ■ — y" — I- )sinh 2 /3 £ Z - 

— I ^ — cos2/3 , Z — -^-sinh2/3 £ Z + ^ cos2/3 £ Z — 


Z 2 


2 t 3 


/3 £ Z cos2/3 £ Z + ( )sinh 2 /3 £ Z + ^ 
2 4pi 5 


( 26 ) 


The equations (25) and (26) can alternatively be derived by replacing /3 £ = 


a,- 


The kinetic energy of the reflector is 
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T 2 = (1/2) m 1 VjV X} -Tn 1 Vja m U )<*+m 2 Vjd_(L >-(1/2) m 2 (± T a_ r (L )5(Z, )<o 

+/n 2 co r 5(L)i.(L)+(l/2) m 2 i 7 '(L)i.(L)+(l/2) ft 7 1 2 Q (27) 

where m 2 is the mass of the reflector and 1 2 Is the mass moment of inertia matrix 
of the reflector. The deflection vector d_(L ) at the mass center of the reflector is 
given as 


d(L) = 


u x {L)- r y u ) 
UyiL) + r x u ^ 
u x tL )r x + u y tL )r y 


(28) 


and the position vector from the point of attachment to the reflector mass center is 
given by 


Thus, 



(29) 


tfU) = 


u x (.L ) — r y zi^(L) 
u y (.L) + r x u ^(L ) 
u x '{L )r x + u y '{L )r y 


(30) 


The angular velocity of the reflector in the inertial co-ordinate system jQ can be 
shown to be 


Cl — 6) + 


U r 


U 


y 

u 


The equation (27) can be simplified as 

T 2 = (1/2) m tVJV* - m.VjaiL )<o+m 2 Yjd_U )+( 1 / 2) m 2 L‘ 


O) j +o> 2 


+m 2 (oJa(.L ~)d_{L )+(l/2) m 2 


Z Z <t>xi U )?< + 

i = U = l 


( 31 ) 
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£ £0 yi (L)<^ (£)?<?, +(l/2)£ r /2£+(l/2)co r /2to (32) 

i=lj = l 

where 

F= 

= £^,'a )«,(0 i^CL^cr) . (33) 

i= 1 i=l i=l 

The kinetic energy of the shuttle, T 0 , is given as 

r 0 = Cl/2)m 1 V 7 ’V.+(l/2) [/i]q. (34) 

where m j is the mass of the shuttle and / 1 is the mass moment of inertia matrix 
of the shuttle. 

The total kinetic energy is given as 

T = T 0 +T x +T 2 (35) 

This can be simplified as 

T = (1/2) m 0 V T V+OL T [//jv+(l/2) <± T j/ 0 ]<*+p£ £ ■ q t 2 +V T a 
+6> r r a+ co^ Q+m 2 Vj d_{L )+m rd_^L )+ 

“'it \ 

m jca T a (Z. )d (Z. )+(l/2) /ft 2 £ 0^(Z. )+0^(Z. ) <ft 2 

[i = l l 

+(l/2)P r /2£+(l/4)p Zp5i?i 2+ I.P6i4i 2 ( 36 ) 

i=l <=1 

where 

m„ = m!+pL +m 2 
H = |pL +m^r_+m 2 ^L}+pLc_ 

1 0 °| 

J 0 = /i+Cl/3) pL 3 0 1 0 +1 2 +J 2 -pLrL-pLrc_-m^rr_-m 7 ra_iL ) 

[0 0 0 
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The term / 2 in this equation can be shown to be: 


Jo = m- 


(r y 2 +L 2 ) ~r x r y 


-r x r. y 

r x L 


r x L 


{r x +L 2 ) r y L 


r y L 


W+r, 2 ) 


The total kinetic energy expression can be further simplified as 


T = Cl/2) m 0 V T V+(jS |h|V+( 1/2) co r 

+w r [A 2 ]i.+Cl/2)i. r |a< 


o>+V r 


where 


= <x_+m 2 d(L ) 


A 2 £= fa+i+m^d )+m ? a (L )rf_(L ) 


A * = 


0 

pL +m 2+p si +p 6 t 
0 


[(*'(z.)) r / 2 (*-a))] 


In this equation 


(37) 


0'CL)) 7 ^ 


0 lx U) 0 0 

0 0 ly 'U) 0 
0 0 <i>\^L ) 


dh/U) 0 0 

0 0 iy tL) 0 
0 0 <t> i ^{L') 


Here i=2,3, n. The number n indicates the total number of flexible modes con- 

sidered. 


Equations of motion 
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Lagrange’s equations of motion for the case of independent generalized co- 
ordinates q k are 


d ST _ JT_ 
dt a q h 8 & 


= Qk ~ 


M. 

8 ?* 


U = 1,2, ,71 ) 


where, T — T{g_gJ is the kinetic energy 
U = U Ca.) is the potential energy, and 

Q k are the generalized forces arising from nonconservative sources. 


(38) 


The generalized co-ordinates are: 

R x ,R y ,R z ~ position of orbiter mass center relative to inertial frame origin. 
0 i, 0 2 >@3 — roll, pitch and yaw angles of orbiter. 
q >?n — modal deformation co-ordinates for the beam. 


The previous kinetic energy expression developed in equation (37) is given in 
terms of nonholonomic velocities V and c^, and generalized velocities q_. Using the 
notation T(V .as .q ) for this kinetic energy expression and T for kinetic energy 
expression in terms of generalized velocities, the equations of motion are developed. 
Thus, equation (37) is rewritten as 

f = (1/2 }m 0 V T V + co 7 [h]v; + (l/2)<o r j/ 0 ]o> + V T (37) 

+ “ r |A 2 |i. + (l/2)i. r |A 3 |i. 

(a) Translational Equations 

From the chain rule applied to equation (37) using equation (8), one gets 


a t 


ar 

bRx 


ev’i 

a^ 

= c T 


9 Ry 


8^2 



af 

a r z 


8^3 


( 39 ) 


Also, the generalized forces are CFO ) where 


£0)= £,0)+F 2 0) 


( 40 ) 


Eo (* ) represents the force applied at the orbiter mass center and F 2 0 ) represents 
the force applied at the reflector mass center. From Lagrange’s equations 


d 

<)T 

+ CC T 

c\T 

dt 

6Z 


az 


= FO) 


(41) 


and from equation (37) 


3L 

BY. 


m 0 V—Hto+A 1 £_ 


Substituting equation (42) in (41), 

m 0 V-H(*+A i ci= -CC T (m 0 V-H(»+A 1J i') + F(f) 


(42) 


(43) 


This can be rewritten as 


m 0 V_ — Ho + A ]£.= Ni + Fit) (44) 

where the nonlinear term N x is given as 

= -CC T (m 0 V-Hu L +A 1 <i) (45) 

= — to(m 0 V—H(±+A l£) 

Here, to = CC T . 

(b) Rotational Equations : 

From equation (4) 

to = M t 6 


Again using the chain rule 
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ar 


a^ r 

[ a« r 

a^i 


a^i 

a^i 

C)T 


f)V T 

|i£.|+ .a a. 7 

a®2 


b0 2 

an] ae 2 

ar 


bv T 

aa r 

a^3 


a®3 

1 903 


It can be shown that 


.6— ^ = y? q b£L i=123 


.9— 7 = o) 7 " Af -1 / = 1 2 3 . 

ae, - ae* 



* TcS §7 



ar| _ 
ai| 

£Tc ^r 

|^zi| + 
l av: ] 

“■'"“asr 


f c *k 




From equation (37), 


and as before 


H =HV+I'0.+ A x i_ 


— m oY. — Ho + A^. 


I 
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Using the Lagrange’s equations 

- MG (51) 

00 

where G is the net moment about the mass center of the orbiter with respect to 
the body-fixed frame. It is given as 

G = Go + (r +a )x.F-> (52) 

Go is the external moment applied about the mass center. Eqation (51) can be 
simpified by substituting equations (42), (49), and (50) together with the relation- 
ship developed in (46) as 


d_ 

dt 


bT 

00 


HV+I 0 v+A £ = G+N 2 


(53) 


where the nonlinear term N_ 2 ls given as 


— A/f-l 


TV 2 = M 


v - ch ^ 

v t c^£— 
- 002 

V T C ^ 


00 , 


3L 

0E 


+M 


-l 


“■'""‘isf 


-Ml 


bT 

0« 


. (54) 


(c) Vibration Equations of the Beam 


Since T in equation (37) is given in terms of g. which is a vector of general- 
ized velocities, 

$L = $L 
bl bL 

and 


bT 

bL 


= A\V+Alto+A 32 . 


(55) 
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The potential energy in the beam is given by 

U = (1/ 2~)gJ Kcl (56) 

where the stiffness matrix K is given as 


K = 


\ o 

hi 

0 \ 


Jr.. 
* It 


and 


(57) 


k u =Eltf 


f (s )ds + f <t>£(s)ds 


<t>li(s ~)ds 


G represents the modulus of rigidity of the beam and j3^ = 
the mass per unit volume (mass density) of the beam. Thus, 


Do>i 2 


r i> 


where D is 


bZ. 


= K<l , 


Using the Lagrangian Equations (38) and assuming that F_i = Q. > 


A\V_ + A 2 & + A —Kg_ . 

(d) Slewing Equations 


(58) 


(59) 


If it is considered to perform a slew maneuver about an arbitrary axis X. and 
the slew angle to be then the slew maneuver can be expressed in terms of four 
Euler parameters. These four Euler parameters are defined as 


6 = 


*2 


= X sin-|- 
— 2 


(60) 
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and their derivatives with respect to time are given as 


(61) 




d € 4 

dr 


1 

CO’ € 

2 -- • 


(62) 

(63) 


If a slew maneuver is considered to be purely rotational, then the transla- 
tional velocity and acceleration can be shown to be negligible during the slew 
maneuver and only the rotational and yibration equations are reqired for the 
analysis and they are simplified by setting V = 0.in both (53) and (59) and are 
written as follows 


I 0 <A+A& = G(r) + AT 2 (co) (64) 

Az&_+ Ag_= —Kg_ . (65) 

Thus equations (62) - (65) completely represent the dynamics of the slew 
maneuver. These equations are nonlinear and coupled including both the rigid- 
body dynamics and the dynamics of the flexible appendage with kinematic non- 
linearities. It is important to note that the nonlinear term /V 2 (<o) is dependent on 
the rotational velocity and as a result determined by the slew maneuver rate. Thus 
the basic slew maneuver stretegy has to be developed before this term can be 
linearized. 

(e) Vibration Equations of the Beam with Damning 

If damping is included in the derivation of vibration equations of the beam, 
then the damping effect can be expressed in terms of frictional forces. These are 
nonconservative, retarding forces and are assumed to be proportional to the gen- 
eralized velocities. In deriving the vibration equations by means of Lagrange’s 
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equations, the following function is introduced 


1 n n 

Zb, 


n is i-* 

Z i=lj = l 




C66) 


It also has a positive definite quadratic form similar to the kinetic and poten- 
tial energy expressions. 


With this definition, Lagrange’s equations assume the form 


d_JtL 6 r b F d_ -q ML (jt = l,2,. — ,n ) (673 

dt 3 ft 3ft 3ft “aft 


Again, as before 


ei. 


= A^V + Alfa + A& 


(55) 


and 


b2- 


= Kg_ 


(58) 


and it can be seen from (66) that 


/ v 



where the damping matrix B is symmetrical and is given as 


( 68 ) 


^11 ^12 • • ^1 n 

5 21 5 2 2 • • 6 2n 


( 69 ) 


1 ^Ji2 
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The vibration equations are given as 


A{V + + Asl + Bl= - Kg. . 

The slewing equations (64) and (65) would be modified as 


( 70 ) 


7 2 <o + A-^_= Git ) + N 2 (<o ,g ) 
Afco + A3l + Bl= - Kg. 


Nonlinear Term in the Rotational Equations 


(71) 

(72) 


The nonlinear term N 2 in the rotational equations (64) and (71) during the 
slewing maneuver is simplified as 


where 


N 2 = M~ l 




- M 


7 0 <o + AM 


(73) 


[o o o 

o y i 1 J 


(74) 


<i) T M — I (— co,sin0 2 cos 2 0 3 +<<> 2 sin0 2 sin0 3 cos0 3 ) (a> j sin0 2 sin0 3 cos0 3 

— Q 9 2 cos 0, l 


— « 2 sin0 2 sin 2 0 3 ) (ca 1 cos0 2 cos0 3 -o) 2 cos02Sin0 3 ) J (75) 


0 i r M -i^M. = _J_ 

— 00 3 COS0 2 


(co 2 cos0 2 ) (— O> 1 COS0 2 ) 


(— Ci> 1 sin0 2 cos0 3 + co 2 sin0 2 sin0 3 + a> 3 cos0^)J . (76) 
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Since the transformation matrix, M , is a function of 0 2 and 0 3 , the time 
derivative of M can be expressed by the chain rule as 




From equation (5) 




(— sin0 2 cos0 3 )0 2 (sin0 2 sin0 3 )0 2 (cos0 2 )0 2 
0 0 0 
0 0 0 




(— cos0 2 sin0 3 )0 3 (— cos0 2 cos0 3 )0 3 0 
(cos0 3 )0 3 (— sin0 3 )0 3 0 

0 0 0 


Substituting these equations (78) and (79) in (77) 


(77) 


(78) 


(79) 


M 


(— sin0 2 cos0 3 )$ 2 +(— cos 0 2 sin0 3 )0 3 (sin0 2 sin0 3 )0 2 + (— cos0 2 cos0 3 )0 3 (cos0 2 )0 2 
(cos0 3 )0 3 (— sin0 3 )0 3 0 

0 0 0 

(80) 


From equation (4), this can also be expressed as 


(sin0 2 cos0 3 )(o> 3 cos0 2 sin0 3 (sin0 2 sin0 3 )(o> iCOS0 2 sin0 3 

+6) 2 cos0 2 cos0 3 )+(— cos0 2 sin0 3 ) + 6> 2 cos0 2 cos0 3 )+ (— COS0 2 cos0 3 ) 
(— o> 1 sin0 2 cos0 3 +to 2 sin0 2 sin0 3 (— <i) 1 sin0 2 cos0 3 +o> 2 sin0 2 sin0 3 

+OJ 3 COS0 2 ) +O> 3 COS0 2 ) 


M = 


1 

COS0 2 


(cos 0 3 )(— <o 1 sin0 2 cos0 3 
+co 2 sin0 2 sin0 3 +o> 3 cos0 2 ) 


(— sin0 3 )(— 6> xsin0 2 cos0 3 
+ 6) 2 sin0 2 sin0 3 + to 3 cos0 2 ) 


0 


0 
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COS0 2 (^1 COS0 2 sin0 3 + 61 2 COS0 2 cos 0 3 ) 


0 

0 


( 81 ) 


Also, M 1 is given as 


M~ l - 


COS0i 


cos0 3 cos0 2 sin0 3 — sin0 2 cos0 3 
— sin0 3 cos0 2 cos0 3 sin0 2 sin0 3 
0 0 cos0 2 


Thus, the nonlinear term A r 2 can be rewritten as 


N 2 = A 3 ( <o.0) |/ 0 61+A 


( 82 ) 


Where the term A 3 is 


A 3( 01 . 0 ) = M 1 


0 



- M 




N 2 — A 3( 01 . 0 ) / ft 6i + A ?( 6i ,0 ) A ->q 


= A/ 6i . 0 ) + A 5 ( 6i , 0 ) j ( 83 ) 

where A 4 depends on the rigid-body slewing and is nonlinear in terms of co and 0 . 
The second term relates the coupling between the rigid-body slewing and the flexi- 


ble modes. 
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4. NUMERICAL DATA 


The analytics developed in the previous section are utilized together with the 
basic SCOLE data [1] and the three dimensional linear vibration analysis [4] to 
generate the following numerical data. 

m l = 6366 A6slugs. ; m 2 — \2A2slugs. ; p= 0.0955 slugs I ft.; L = 130/f. 
G+ = 7.2E+Slb I ft 2 ; (£7 ) x = (El ) y = (El) = AE+Vb-ft 2 ; 


0.036 
r_ = -0.036 
-0.379 


c_= 0 

-65.0 

905443.0 0.0 145393.0 

I I = 0.0 6789100.0 0.0 

145393.0 0.0 7086601.0 

18000.0 -7570.0 0.0 

1 2 = -7570.0 27407.0 0.0 

0 0.0 27407.0 

The three dimensional vibration analysis is given in terms of the first ten 
modal frequencies and mode shapes in table 1. Here, 

OtiS &4S (XjS OtjS 

<l> xi is ) = A xi sin— — h B xi cos— — VC xi sinh— — h D xi cosh— — 

i j i ■/ j 


VCi& OliS OCjS 

<t>yi (s ) = A yi sin— — I -B yi cos— — \-Cyi sinh— — \-D yi cosh 

JLt J—* J-* 


&\pi ) A yjji sinoiyjj} j B cosof ^ 
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i_ 

4 


dl 1 ^ It 

G 

Using these data the following matrices are obtained. 

1216640 15167.53 - 115118.9 
I 0 = 15108.34 7083005 - 52474.84 

[-115096 - 52503.9 7131493 
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THREE DIMENSIONAL MODE SHAPE CHARACTERISTICS 

MODE No. 

1 

2 

FREQ. (Hz.) 

0.27804240E+00 

0.31357296E+00 


MODE No. 


FREO. (Hz.) 


MODE No. 


FREO. (Hz.) 


0.12012084E+01 
0.16282665E+00 
-0.19670286E+00 
-0.16983450E+00 
0.196 16259E+00 
-0.10274618E-01 
0.57579 133E-02 
0.11810057E-01 
-0.57220462E-02 
0.19360955E-01 
-0.50748354E-01 
0.1 397801 8E-04 


3 


0.8 1300 189E+00 


0.20540387E+01 
0.40868 188E-01 
-0.6 195884 5E-01 
-0.41309992E-01 
0.61880796E-01 
-0.22438404E-01 
0.36509234E-01 
0.2439044 7E-01 
-0.36464 758E-01 
0.566 11842E-01 
0.9269890 IE-01 
-0.87320799E-05 


5 


0.20536300E+01 


0.32645546E+01 
0.99278129E-01 
-0.92344553E-01 
-0.99442 145E-01 
0.9222580 IE-01 
-0.57396019E-01 
0.53976008E-01 
0.581 14853E-01 
-0.53906980E-01 
0.1 430006 2E+00 
-0.16588614E-02 
0.61861804E-07 


0.1275651 8E+01 
0.3885529 IE-02 
-0. 1499838 7E-01 
-0.4332101 8E-02 
0.14985820E-01 
0.142 1978 1E+00 
-0.22695797E+00 
-0.19283 105E+00 
0.2264456 1E+00 
0.21835058E-01 
0.311 15282E-01 
-0.7599233 7E-05 


0.11856099E+01 


0.24804687E+01 
0.80641 794E-01 
-0.67233377E-01 
-0.809 13938E-01 
0.67106316E-01 
0.13728679E+00 
-0.1 174693 2E+00 
-0. 14085209E+00 
0.1 172505 7E+00 
0.82557693E-01 
-0.16158934E-03 
0. 1043771 8E-07 


6 


0.497 16090E+01 


0.497 16090E+01 
0.45739784E-01 
-0.4636558 IE-01 
-0.45763 106E-01 
0.463 29676E-01 
0.786 12940E-01 
-0.7995285 3E-01 
-0.789 14485E-01 
0.79891039E-01 
0.33165303E+00 
-0.93394833E-05 
0.1501721 IE-09 
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THREE DIMENSIONAL MODE SHAPE CHARACTERISTICS 


MODE No. 

7 

8 

FREQ. (Hz.) 

0.55157833E+01 

0.1228 124 9E+02 

a 

0.5350 1560E+01 

0.79833305E+01 


0.8131 1804E-01 

0.4483506 IE-01 

B x 

-0.82056569E-01 

-0.4483491 4E-01 

c x 

-0.81344923E-01 

-0.44840508E-01 

D x 

0.81997259E-01 

0.44813000E-01 

Ay 

-0.47 145439E-01 

0.77404756E-01 

By 

0.477035 90E-01 

-0.77465629E-01 

Cy 

0. 47289807E-01 

-0.77475327E-01 

Dy 

-0.47669155E-01 

0.77427782E-01 

OiyU 

0.384081 10E+00 

0.85518143E+00 

A* 

-0.23855560E-02 

0. 1583037 IE-05 

B & 

0.33122041E-07 

-0.9871501 7E-11 

MODE No. 

9 

10 

FREQ. (Hz.) 

0.12890442E+02 

0.23679520E+02 

a 

0.81789349E+01 

0.11 08534 7E+02 

A x 

0.78743585E-01 

0.44348498E-01 

B x 

-0.78755259E-01 

-0.44367373E-01 

c x 

-0.78752483E-01 

-0.4435051 IE-01 

D x 

0.7871 7693E-01 

0.4435 1763E-01 

Ay 

-0.4556 9244E-01 

0.76707490E-01 

By 

0.45609474E-01 

-0.76762782E-01 

Cy 

0.45607884E-01 

-0.76733612E-01 

Dy 

-0.4558 7726E-01 

0.76735779E-01 

a * 

0.897601 45E+00 

0.16488784E+01 

A* 

0.94995483E-03 

-0.51 10595 7E-06 

B 

-0.5643 7766E-08 

0.16528495E-11 















— 0 . 89042 E — 1 
0 . 11263£0 
0 . 32737£+2 
— 0 . 6392£— 1 
- 0 . 14526£0 
0 . 7925 E — 2 
— 0 . 8369 E — 1 
- 0 . 2829£-2 
0 . 89767 E -1 
0 . 1866 E — 2 


- 0 . 14067£0 
— 0 . 1471£0 
— 0 . 6392 £ — 1 
0 . 2547£+3 
0 . 1908£0 
— 0 . 4278 £— 1 
- 0 . 761 15 £— 1 
0 . 1543 £— 1 
0 . 2859 E — 1 
— 0 . 9067£— 2 


0 . 1457£0 
— 0 . 5518£— 1 
— 0 . 14526£0 
0 . 1908£0 
0 . 8103£+3 
— 0 . 2570£— 1 
— 0 . 12912£0 
0 . 9222£— 2 
0 . 461 1£-1 
- 0 . 594 7£ -2 


0 . 45879£+2 
0 . 36305£— 1 
— 0 . 89042£— 1 
— 0 . 14067£0 
— 0 . 1457£0 
0 . 1914£— 1 
0 . 84597£— 1 
— 0 . 6893£— 2 
— 0 . 426 9 £—l 
0 . 4204£— 2 


0 . 36305£— 1 
0 . 621 1£ +2 
0 . 11263£0 
— 0 . 1471£0 
— 0 . 551 8 £—l 
0 . 19839£— 1 
0 . 3935£— 2 
— 0 . 7165£— 2 
0 . 5969£— 2 
0 . 41227£— 2 


— 0 . 6893£— 2 
— 0 . 7165£— 2 
— 0 . 2829£— 2 
0 . 154 3 £— 1 
0 . 9222£— 2 
— 0 . 2089£— 2 
— 0 . 37286£— 2 
0 . 1342962£+8 
0 . 142 1£— 2 
— 0 . 4427£— 3 


— 0 . 4269£— 1 
0 . 596 9£— 2 
0 . 89767£— 1 
0 . 285 9 £-l 
0 . 461 1£-1 
— 0 . 3955£— 2 
— 0 . 385 9 £—l 
0 . 142 1£— 2 
0 . 2095672£+8 
— 0 . 9108£— 3 


0 . 4204£— 2 
0 . 4127£— 2 
0 . 1866£— 2 
— 0 . 9067£— 2 
— 0 . 5947£— 2 
0 . 1227£— 2 
0 . 2397£— 2 
— 0 . 4427£— 3 
— 0 . 9108£— 3 
0 . 8662547 £+ l 0 


0 . 1914£— 1 
0 . 19839£— 1 
0 . 7925£— 2 
- 0 . 4278£— 1 
- 0 . 2570£— 1 
0 . 23209£+5 
0 . 10383£— 1 
- 0 . 2089£— 2 
— 0 . 3955£— 2 
0 . 1227£— 2 


0 . 845 97 £—l 
0 . 3 93 5£— 2 
— 0 . 8369 £— 1 
— 0 . 76115£— 1 
— 0 . 12912£0 
0 . 10383£— 1 
0 . 5556 l £+5 
— 0 . 37286£— 2 
— 0 . 3859£— 1 
0 . 239 7£— 2 


A\ 


— 0 . 2133821£0 
0 . 3808921£+3 
— 0 . 1808478£+3 
0 . 1423380£+3 
— 0 . 2416743£+2 
— 0 . 6802273£0 
0 . 2784792£+2 
0 . 78428 18£+1 
— 0 . 2694455£ +2 
— 0 . 9225328 £— 1 


— 0 . 3687057£+3 
— 0 . 3030935£+2 
- 0 . 13185 96£+3 
— 0 . 1135851 £ + 1 
0 . 574383£+2 
0 . 3104929£2 
0 . 6651585£+2 
— 0 . 1930097£+2 
— 0 . 5544252£ +2 
0 . 1594045£+2 


— 0 . 7253901£— 1 
— 0 . 8427658£— 1 
— 0 . 125799£0 
— 0 . 2367351£— 1 
— 0 . 9150328£— 1 
- 0 . 384 3062 £-l 
0 . 596075£— 1 
— 0 . 4363533£— 2 
— 0 . 4200623 £— 1 
— 0 . 1626004£— 1 
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The stiffness matrix K is calculated using equation (57) and the mode shape 
coefficients given in Table 1. This matrix is a diagonal matrix and is represented in 
terms of the diagonal elements as 


K = 


* u = 0.28202 17£0 
* 2;2 = 0.35746 92£0 
*3 3 = 0.2412807£1 
* 4 , 4 = 0.52851 16£1 
* 5 , 5 = 0.1588654£2 
* 6 , 6 = 0.8573860E2 
* 7 ,7 = 0.1 1461182s 3 
* 8 , 8 = 0.5686101E3 
* 9 , 9 = 0.62545982s 3 
*10,10= 0.21 14612£4 


The damping matrix B used for this analysis is a diagonal matrix and for 
damping ratio £ = 0.003, it is calculated to be 


B = 


b x> i = 0.96859642s— 3 

b 2 ,2 = 0.10886082s— 2 

b 3i3 = 0.28340 16£— 2 
b 4 , 4 = 0.4256808£— 2 
b 5 , 5 = 0.7387177£-2 
b 6>6 = 0.17190142s — 1 
b 7 , 7 = 0.1984237£-1 
2 > 8 , 8 = 0.442 1234£-1 
b 9> 9 = 0.46334342s— 1 
* 10 , 10 = 0.852764 7£ — l 
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APPENDIX 

The following is a summary of transformations between inertial frame and 
body-fixed frame. Here, s* and q (z = 1 , 2 , 3 ) denote sinQj and cos 0 j (z = 1 , 2 , 3 ) 
respectively. 


(a) Space-three Angles 


C = 


C 2 C 3 c 2^3 ~ S 2 

SiS 2 C 3 —‘ S 3 C X S x S 2 S 3 +C 3 Ci s 1 c 2 

c \S 2 c 3+s 3S 1 c^^a— c 3 s x c x c 2 


M t 


1 0 -s 2 

0 C ! S X C 2 
0 —Si c x c 2 


(b) Space-two Angles 


C = 


C 2 S 2 S 3 ~ S 2 C 3 

s x s 2 —s 1 c 2 s 3 +c 3 c l 5^203+530! 

C X S 2 -C X C 2 S 3 -C 3 S 1 C jC 2 C 3 S 3 S 1 


m t 


10 c 2 

0 cj s x s 2 
0 —Si c is 2 


(c) Body-two Angles 


C = 


c 2 s x s 2 -c x s 2 

S 2 S 3 -S X C 2 S 3 +C 3 Ci C 1 C 2 S 3 +C 3 S 1 
S 2 C 3 -S 1 C 2 C 3 -S 3 C 1 CiC 2 C 3 -S 3 Si 


m t 


c 2 0 1 

S 2 S 3 C 3 0 

S 2 C 3 —S 3 0 





~2 *B2 *B1 


Z B3" Z B2 


(c) 


(d) 


(a) Axes in reference position (b) First rotation-about x axis 
(c) Second rotation-about y axis (d) Final rotation-about z axis 


FIGURE 2 
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Figure 3- Position Vectors in Inertial Frame 



Figure 4- Vectors in Body-fixed Frame 
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16. Abstract 

This is the first report of a set of two reports on the dynamics and control of 
slewing maneuvers of NASA Spacecraft Control Laboratory Experiment (SCOLE) article. 
In this report, the dynamics of slewing maneuvers of SCOLE are developed in terms 
of an arbitrary maneuver about any given axis. The set of dynamical equations 
incorporate rigid-body slew maneuver and three-dimensional vibrations of the 
complete assembly comprising the rigid shuttle, the flexible beam, and the reflecto 
with an offset mass. The analysis also includes kinematic nonlinearities of the 
entire assembly during the maneuver and the dynamics of the interaction between 
the rigid shuttle and the flexible appendage. The final set of dynamical equations 
obtained for slewing maneuvers are highly nonlinear and coupled in terms of the 
flexible modes and the rigid-body modes. 


The equations are further simplified and evaluated numerically to include the first 
ten flexible modes and the SCOLE data to yield a model for designing control 
systems to perform slew maneuvers. 












